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Objective 

This  work  derived  the  mathematics  required  in  the  RSCEM  (Radar  Simulation 
Concept  Evaluation  Model)  model  for  closed-form  answers  for  the  probability  of 
detection  and  false  alarm  for  Swerling  class  targets  in  clutter  and  noise  and  targets  in 
compound-K  clutter  plus  noise.  In  all  cases,  the  results  are  applicable  to  the  N-pulse 
coherent  detection  problem.  For  the  case  of  a  target  in  “spikey”  clutter  plus  the  noise  the 
advantage  of  multiple  pulse  integration  is  significant. 

Funding  was  provided  under  contract  #N66001-94-D-0009. 
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1.  Introduction 


Target  detection  in  Gaussian  noise  is  an  old  problem  and  there  too  many 
references  to  cite  them  all  here.  Similarly,  there  are  several  published  results  for  target 
detection  in  Gaussian  noise  and  log-normal  clutter  (Trunk  and  George,  1970,  Schleher, 
1975).  Although  some  of  the  analysis  in  this  report  is  derived  elsewhere,  we  have  tried  to 
present  a  unified  approach  to  the  problem  of  target  detection  in  clutter  plus  noise  to 
include  all  the  standard  Swerling  target  cases  and  the  relatively  new  compound-K  clutter 
probability  density  function,  as  discussed  by  Ward  et  al.(1990).  In  the  following  analysis, 
many  of  the  derivations  are  new  and  yield  closed  form  results  for  the  N-pulse  case  for  the 
probability  of  detection,  and  probability  of  false  alarm  for  Swerling  class  targets  in 
Gaussian  noise,  assuming  square-law  detection.  All  the  derivations  are  not  new; 
however,  to  the  authors  knowledge,  they  do  not  appear  in  a  single  reference  for  all 
combinations  of  Swerling  targets  and  clutter  statistics,  including  compound-k.  The  report 
considers  the  following  theoretical  problems: 

1)  The  use  of  the  characteristic  function  to  treat  the  problem  of  N-multiple  pulse 
integration,  for  SC-0,  SC-1,  SC-2,  SC-3,  SC-4  and  compound-k  clutter  cases. 

2)  The  use  of  a  basic  rule  in  probability  theory  for  removing  “left”  and  “right” 
variables  in  conditional  pdfs  to  treat  the  problem  of  target  plus  clutter  in  noise,  and  then 
target  in  clutter,  in  place  of  multiple  convolutions  of  the  individual  pdfs  representing  the 
sum  of  the  random  variables  target,  clutter  and  noise. 

3)  The  use  of  the  recently  reported  compound-K  pdf  for  treating  “spikey”  clutter 
arising  in  wide-band  radar  applications. 

4)  FORTRAN  77  algorithms  for  all  Swerling  cases  for  target  plus  noise,  and 
Swerling  case  #1  and  #2  targets  in  compound-k  clutter. 
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2.  Probability  Density  Functions  for  a  Target  in  Clutter  and  Noise 
2.1  General  Case 

A  detection  decision  is  to  be  made  during  each  N-pulse  beam-dwell  interval.  In 
the  general  case,  the  target  and  the  clutter  can  vary  from  pulse-to-pulse  and  from  scan-to- 
scan.  An  extension  of  the  usual  two-components  representing  the  complex  phasors  for 
the  I  and  Q  output  of  the  received  signal  to  include  clutter,  after  post  detection  is  the 
following  three  component  model 


Re'®  =  T  e'“  -H  -h  re'^ 


(1) 


where  t  e'“  is  the  complex  phasor  representing  the  target,  is  the  complex  phasor 


representing  the  clutter,  and  re  is  the  complex  phasor  representing  thermal  noise.  A 
vector  diagram  for  the  three  phasors  is  shown  in  figure  1,  where  the  detection  is  made  on 
the  basis  of  Re  in  equation  (1).  Note,  the  symbol  CO  in  equation  (1)  should  not  be 
confused  with  the  same  parameter  when  used  in  the  characteristic  function. 


Q 


Figure  1.  Phasor  diagram  for  three  components  of  the  received 
signal. 


If  the  target  is  immersed  in  Gaussian  noise  and  clutter,  and  this  target  is  to  be  detected 
using  an  N-pulse  radar  burst,  the  return  can  in  equation  (1)  can  be  written  as 


Re*  =  pe”  +  re" . 


(2) 


where  the  target  plus  clutter  is  decomposed  as 
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(3) 


The  pdf  for  the  return.  Re'®,  can  be  written,  using  the  basic  rule;  to  remove  a  number  of 
“left’  variables,  we  integrate  with  respect  to  them.  To  remove  a  number  of  “right” 
variables,  we  multiply  by  their  conditional  density  with  respect  to  the  remaining  “right” 
variables  and  integrate  (Papoulis,  1965,  p.236).  For  example,  for  two  “right”  removals 
we  have. 


1^4  )  =  JJ  kz  ’  ^3  >  ^4  )p{^2  ’^3^4  )dx^dx^  .  (4) 


For  the  case  of  a  target  in  noise  alone  (no  clutter  in  equation  (1) )  standard  Swerling  cases 
use  a  Rayleigh  pdf  for  the  noise  component  given  by 


(5) 


The  distribution  of  R  in  equation(l)  with  no  clutter;  i.e.,  where  =  T  is 
Rice  given  by 


a: 


+p^)/2aj 


o? 


(6) 


Assume  a  square-law  detector.  The  actual  detector  type;  i.e.,  linear,  half-wave,  square- 
law,  full-wave  square-law  introduces  only  minor  differences  in  the  correlation  function  of 
the  received  signal,  i.e.,  all  correlation  outputs  are  proportional  to  the  mean  square 
received  power  with  second  order  differences  depending  on  the  time  shift.  Make  the 
following  substitutions: 
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RdR 


R^  n 


(7) 


and  temporarily  treat  p  as  though  it  were  a  constant,  and  R  is  the  detected  amplitude. 

Substituting  equation  (7)  into  equation  (6)  gives,  for  each  pulse,  in  the  SC-0,  SC-1,  SC-2, 
SC-3,  and  SC-4  cases. 


0  ,y<o’ 


(8) 


where  the  pdf  in  equation  (8)  is  conditioned  on  x.  DiFranco  and  Rubin  (p.370,  1980) 
show  that  there  is  negligible  difference  in  detection  performance  (as  measured  in  terms  of 
the  probability  density  function)  between  a  square-law  envelope  detector  and  a  linear 
detector  (on  the  order  of  0.1  dB).  Blake  (p.44,1986)  also  states  the  difference  between 
square-law  and  linear  detection  in  terms  of  the  required  signal-to-noise  ratio  for  detection 
is  on  the  order  of  0.2dB  at  most.  As  a  point  of  interest,  the  result  in  equation  (8)  can  be 
shown  to  be  equivalent  to  the  result  in  Nathanson.* 


'  second  edition  (p.  168,  Table  5.6,  entry  (5),  Rice  (power),  1991),  by  making  the 
following  substitutions  in  equation  (8) 

y  =  {\  +  a^)x 
p  =  a4x 

R^I2  pV2 
y  =  — ^ ,  X  =  — 


where  is  the  ratio  of  th^arge  reflection  (target)  to  the  sum  of  small  reflections 
(clutter),  X  is  the  res,  and  x  =  2(7^  is  the  mean  res.  We  note  a  “typo”  in  Nathanson, 

where  the  density  function  for  the  res,  X,  should  read  in  the  argument  of  the  Bessel 
function  as 


J 


o 


=  ^.(2a/^) 
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In  order  to  treat  the  case  of  multiple  pulse  post-detection  integration  gain 
efficiency,  it  is  convenient  to  introduce  the  characteristic  function  associated  with  ^(>^1 
in  equation  (8)  as 


^(co|x)  =  j  p(y\x)e‘'^-dy. 


(9) 


and  substituting  equation  (8)  into  equation  (9)  gives 


4>(0)|4  =  e-']e-*'-‘‘X{24^)dy 

0 


X 

,(1-/0)) 


-e 


(1  -  /co) 


(10) 


1-/(0 


We  will  now  consider  the  N-pulse  detection  case.  For  N-pulse  detection,  we 
detect  the  sum  of  N  different  values,  y, , . . . ,  of  the  received  power  and  noise  as 


+...+ 


f  Rl/2^ 


V 


J 


2o^  2g^ 

n  n 

Exp{Y\x}  =  N^-^  +  N  =  Nx  +  N  =  N(x  +  \) 


(11) 
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where  the  last  line  in  equation  (11)  will  be  obtained  by  an  alternative  method  shortly,  and 
X  is  defined  in  equation  (7). 

If  the  pulse-to-pulse  noise  samples  are  statistically  independent,  then  the 
characteristic  function  of  the  sum  of  the  N  samples  is  the  produce  of  their  individual 
characteristic  functions.  Therefore,  the  characteristic  function  associated  with  Y  is 


(l-mf 


(12) 


A  slight  modification  of  the  order  of  first  averaging  over  p(x)  and  next  taking  the 
product  of  the  individual  characteristic  functions  is  required  for  SC-4.  In  this  case  we 
need  to  average  with  respect  to  one-dominant-scatter  plus  Rayleigh  noise,  before  finding 
the  sum  of  the  N  samples.  The  advantage  of  performing  the  multiple  pulse  analysis  in  the 
“frequency”  domain  as  in  equation  (12)  as  opposed  to  the  probability  domain  is  that  the 
expression  corresponding  to  (12)  in  the  pdf  domain  is  an  N-fold  convolution  of  each  of 
the  N  pdfs  for  each  of  the  N  variables. 

Now,  because  the  target  plus  clutter  fluctuates  from  scan-to-scan,  for  all  Swerling 

1 

class  targets,  we  must  average  equation  (10)  over  the  variable  x  =  -^—.  That  is. 


OO 

— OO 

OO 

=  j[0(co|x)]'';7(Y)fi6c  ,  (13) 

■  OO 

PMdx 

(1  — 1(0)  — 


where  p(^x)  is  the  probability  density  function  for  the  random  variable  X.  In  the  case  of 
Swerling  class  targets,  p{x)  corresponds  to  the  density  function  of  the  target  alone, 
appropriate  to  either  pulse-to-pulse  variation  or  scan-to-scan  variation.  For  the  case  of 
target  plus  clutter  plus  noise,  a  similar  analysis  will  be  used  as  shown  in  section  2.7.  In 
the  case  of  Swerling  class  targets,  once  ^^(Cd)  is  known,  the  probability  density  of 
Y  =  y^+y^-\ - 1-  ,  the  desired  result,  is 


II 


(14) 


271  — 


Equations  (13)  and  (14)  are  extremely  useful,  and  form  a  “key  pair”  and  will  be 
used  in  several  of  the  following  derivations. 

Before  leaving  the  subject  of  the  characteristic  function  associated  with  a  pdf  it  is 
instructive  to  consider  the  meaning  of  the  parameter  CO^  in  the  definition  for  the 
characteristic  function.  First,  this  parameter  is  a  dummy  variable  of  integration.  Second, 
it  is  not  really  frequency,  but  simply  a  variable  used  to  analytically  continue  the 
characteristic  function  into  the  complex  plane  in  order  to  consider  the  moments  of  the 
random  variable.  The  derivatives  of  the  characteristic  function  of  a  random  variable,  Y, 
are  related  to  its  moments  as 


Exp{Y-\x} 


1  d"^ 


r  ddi" 


(0=0 


(15) 


and  the  Taylor’s  series  expansion  (a  Laurent  series  is  a  more  general  analytic  continuation 
for  a  function  of  a  complex  variable)  for  the  characteristic  function  is 


/  ,  .  ,  (i(Oy 

<I>(co|jc)  =  1  + /com,  -CO  m^  +•••  +  — —m^  + 


(16) 


For  the  characteristic  function  given  in  equation  (16)  we  find  for  the  first  two  derivatives 


^  See  remark  on  page  6. 
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_  j  e 


1-/(0 


Jco  0  (1  -  /co)  L 
d<b 


iN  +  iNx 


Ngxk 

1-/0 


d(x) 


(0=0 


= /A^J  (1  +  j£:)/7(x)(/x: 

0 

=  -  A/J  xp(x)dx  +  iN]  {iN  +  iNx)  p{x)dx 


(17) 


(0=0 


and  from  the  second  line  in  equation  (17)  we  see 


Exp{Y\x}  =  A/(x  + 1)  , 


(18) 


which  is  the  same  statement  as  equation  (1 1)  as  of  course  it  must  be.  It  is  also  easy  to 
show  that 


=  -N[1  +  2(x  +  n  +  A/x)] 

cl=^Exp{Y^}-{Exp{Y}y  (19) 

=  3N 
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Before  considering  the  various  Swerling  cases,  it  is  instructive  to  draw  a  picture  of 
target  echo’s  varying  on  either  a  pulse-to-pulse  or  a  scan-to-scan  basis.  Hopefully,  such  a 
picture  helps  explain  the  choice  for  a  particular  pdf  used  in  the  following  sections.  Figure 
2  shows  a  time  series  (or  equivalently,  a  range  profile)  of  received  echo  fluctuations 
under  different  circumstances. 


Scan#l 


Scan  #2 


SC-0:  Steady  target. 


Scan  #1  Scan  #2 

SC-1  and  SC-3:  Steady  target  pulse-to-pulse. 


SC-2  and  SC-4:  Pulse-to-pulse  variation.  For  SC-4,  there  is  a 
dominant  pulse  in  each  scan. 

Fig.2  Various  Swerling  class  time  series  representations  for  a  fluctuating  target. 
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2.2  Swerling  Case  #0  (non-fluctuating  target): 


In  this  case,  we  see  from  figure  2,  that  the  target  is  constant  pulse-to-pulse  and 
scan-to-scan,  and  the  total  received  signal  can  be  written  as 


Re" 


AV2 


re 


i(p 


(20) 


where 


T  e 


tco 


(21) 


and  where  the  average  target  power  to  noise  power  is  the  first  term  on  the  right-hand-side 
of  equation  (20),  and  the  phase  is  constant.  This  case  was  originally  discussed  by 
Marcum  (1960).  The  pdf  for  the  target  in  equation  (20)  is 


p{x)  =  S 


(  A^I2\ 

X - r- 


V 


J 


(22) 


and  substituting  equation  (22)  into  equation  (13)  gives 


1 

(l-icaf 


(23) 


in  agreement  with  Fehlner  (1962,  his  equation  A-4  on  page  24).  Substituting  equation 
(23)  into  equation  (14)  of  our  “key  pair”  gives 
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(24) 


pAy)=^]  ‘ 


AV2 


We  need  to  perform  a  bit  of  algebra  to  write  the  integral  in  equation  (24)  in 
‘closed  form,”  a  goal  throughout  this  report. 


1  -  /CO  =  iu,  /co  =  1  —  iu 


du  =  -dG) 


Ps(.Y)  = 


f  ^ — -du 

J 


27C  -/-oo  (iw) 

^  =  iu,  =  idu 

K) 


-iY+m 


(25) 


The  last  line  in  equation  (25)  is  can  be  written  in  terms  of  a  Laplace  transform  as 


1 


m 


_  ^-(Y+m 


N-] 

Y 


h.X24m) 


(26) 


where  the  last  line  in  equation  (26)  follows  from  transform  pair  29.3.81,  page  1026, 
AMS55,  Handbook  of  Mathematical  Functions,  National  Bureau  of  Standards,  1964. 
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The  probability  of  detection,  for  the  N-pulse  case,  is  then  the  simple  numerical 
integration  of  equation  (26)  as 


]pAY)dY 


(— 

KNb) 


N-\ 
\  2 


(24JM)dY 


(27) 


The  probability  of  false  alarm  is  found  by  setting  ^  =  0  (N.B.,  this  corresponds 
to  setting  x,  the  ratio  of  mean  target  power  to  noise  variance  to  zero),  and  this  is  most 
easily  accomplished  in  equation  (25)  since  (26)  is  indeterminate.  From  (25)  with  ^  =  0 
we  find 


p-Y  1+,- 

ZTLl  \-ioo  ^ 

{N-l)\ 


(28) 


and  the  probability  of  false  alarm  becomes 


]p„(Y,b=mY 


1 

(V-l)! 


(29) 


which  is  easily  evaluated  by  integration  by  parts. 
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ForaP„  =10-'  (y  =8xlog,10sl8.42068),andN=l, 


P„=l-e-‘Je-7.(2Vw)dy, 

0 


and  for  N=2, 


I 

P,  =  1  -  j  /.  {2-j2bY)dY. 

and 


=e-Hi+y), 


and,  for  illustrative  purposes,  figure  3  shows  a  plot  of  Jr  versus - ; —  in 


Fig.  3  Probability  of  detection  versus  signal-to-noise  ratio  for 
Swerling  Case  #0. 


2.3  Swerling  Case  #1: 


For  both  Swerling  cases  SC-1  and  SC-3  we  see  from  figure  2  that 


Re'®  =p  +  re“^  , 


(30) 


where  the  first  term  on  the  right-hand-side  of  equation  (30)  fluctuates  from  scan-to-scan 
but  is  constant  from  pulse-to-pulse.  The  scan-to-scan  Rayleigh  pdf  for  the  normalized 
mean  target  power  for  this  case  is 


p{x')  =  e~^  ,x'>0. 


(31) 


or,  in  terms  of  the  unnormalized  mean  power  to  noise  variance,  , 


pV2  pV2 


X 

a  n 


a 


X 

or; 

dx'  =  —dx 

X 


(32) 


p{x)  =  p(x')^  =  - 


dx 


and  in  going  from  X  to  x'  there  are  two  normalization’s,  first  by  X^,  and  then  the 
normalization  used  in  the  definition  of  X  itself. 

a 

Before  proceeding,  we  show  that  equation  (32)  is  indeed  a  Rayleigh  pdf,  for  the 
target  voltage,  p .  From  the  definition  of  X, 
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,  p  . 

=  — —,dx  =  —d^ 

a  a 

n  n 

p{p)  =  pix)^  (33) 

dp 


Substituting  p{x)  from  the  last  line  of  equation  (32)  into  equation  (13)  gives 


1 

(l-icof 


iN(0 


l-t(0 


^dx 


_ 1 _ 

(l-tCD)'^'’[l-tCD(l  +  A(xJ] 


(34) 


Substituting  equation  (34)  into  equation  (14)  yields  the  probability  density  for  the  N-pulse 
case  as 


1 


1 


J; 


-imY 


CO  = 


2n  (l  +  Nx^ )  i(co„  -  z0)(l  -  mY' 
1 


day 


(35) 


l  +  Nx 


We  now  treat  various  values  for  N. 
case  N=l: 
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J _ 1_ 

2n  (l  + 


d(0 


\+Xa 


l  +  x.. 


and  the  corresponding  probability  of  detection  is 


(36) 


P^  =  lp,(Y)dy  =  e'^‘-, 


(37) 


and  Pf.^  =  e~^-  (from  equation  (36)  with  =  0),  and  both  36  and  37  agree  exactly  with 
Fehlner  (his  equation  15,  page  18). 

N=2: 

Expand  the  integrand  in  equation  (35)  using  a  partial  fraction  expansion  as 


1 


B 


(co^  -  l(o)(l  -  ico)  co„  -  tco  1  -  /(O 
A  -  fcoA  +  5(D„  -  fco5  =  1 
A  =  -B 
1 


(38) 


A  = 


1-co,. 


and  substituting  equation  (38)  into  (35)  gives 


a(5')  = 


1  1 
(l-(0.)(l  +  2x.) 


(39) 


and  the  corresponding  probability  of  detection  is 
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(40) 


14 


r„ 

l+2x. 


where  the  last  line  in  equation  (40)  is  valid  for  signal-to-noise  ratios  greater  than  10  (10 
dB).  From  the  mathematical  principle  of  induction  it  is  readily  demonstrated  that  for 

arbitrary  N, 


1 


N  W-I 


ik 


e 


1+A!x„ 


(41) 


a  very  simple  closed  form  result.  The  corresponding  probability  of  false  alarm  is 
obtained  from  equation  (41)  by  making  the  substitution  =  0,  obtaining 


(42) 


Figure  4  shows  a  plot  of  versus  X^ 


dB  for  iV  =  1,10. 
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2.4  Swerling  Case  #2: 


From  figure  2,  we  see  the  received  signal  is  given  by 


Re'®  =  pe''"  +  • 


(43) 


where  the  first  term  on  the  right-hand-side  of  equation  (43)  fluctuates  from  pulse-to- 
pulse  and  scan-to-scan.  In  this  case  the  target  is  Rayleigh  (the  form  given  by  Nathanson 
is  exponential  for  the  power  which  is  Rayleigh  for  the  voltage)  and  immersed  m  Rayleigh 
noise  and  the  joint  probability  density  of  T  and  ({)  is  given  by 


r 

24o’+a^) 


(44) 


where  is  the  target  variance  and  is  the  noise  variance.  Since  we  are  again 
interested  in  obtaining  the  probabiUty  density  of  the  square  of  the  envelope  in  equation 
(43)  it  is  convenient  to  define 


/  = 


RV2 


1-+ 


(Jr 


RV2 


(45) 


pV2 


X 


a 


1  +  JC . 


where  equation  (45)  uses  the  same  form  for  the  definition  of  x'  as  in  equation  (7)  for  the 
normalized  total  power.  The  renormalization  of  X  in  (45)  corresponds  to  a  simple 
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substitution  in  all  the  derivations  which  follow.  Equation  (45)  implies  that  the  target 
variance  plus  noise  variance  are  embedded  (or  included)  in  the  pdf  in  equation  (8)  and  do 
not  need  to  be  included  again  in  the  pdf  for  the  target  random  variable,  X.  Thus  the  pdf  is 

a  delta  function  given  by 


p(x)  =  5(x). 


(46) 


The  normalization  in  the  first  line  in  equation  (45)  is  motivated  by  considering  the  total 
received  power  in  each  pulse.  From  equation  (43),  i.e.. 


Re'®  =  pg''"  +  re'* 
lR|^=lpe''^  +  re'*| 

=  +  r^  +  2pr  cos(\}/  -  <|)) 

Exp{R^ }  =  £jcp{p' }  +  Exp{r^ }  +  2£xp{pr  cos(\jf  -  (j))} 


(47) 


For  N-pulse  detection  of  a  Swerling  case  2  target,  where  the  target  is  varying  from 
pulse-to-pulse,  equation  (13)  is  replaced  by 


r'  =  y;+y.'+-"+3^^ 

and,  in  the  probability  of  detection  is  replaced  with  Y^  .  Substituting  equation  (46) 
into  equation  (13)  yields  the  characteristic  function  for  the  SC-2  as 


- L— 

(1  —  ico)  0 
1 


(49) 
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in  agreement  with  DiFranco  and  Rubin  (Equation  1 1.3- 17a,  page  404).  We  will  see 
shortly,  that  SC-2  and  SC-4  differ  only  in  p{x).  Substituting  (49)  into  equation  (14) 
gives  the  N-pulse  probability  for  Swerling  case  #2  as 


(iV-1)!  ’ 


(50) 


in  agreement  with  Fehlner  (his  equation  17,  page  19). 

The  probability  of  detecting  a  Swerling-2  target  in  noise  is  the  probability  that 


r  > 


Y 

_ o 

l  +  x 

a 


(51) 


or 


(52) 


and  subsfltuting  equation  (50)  into  equation  (52)  gives  the  foUowing  probabiUty  of 
detecting  a  Swerling  #2  target  as 


r 


Po=e 


l+El.  N-\ 

"■s 


m-0 


\ 


14 


J_ 

m! 


(53) 


another  closed  form  result.  The  probability  of  false  alarm  is 


-y„ 

e  . 
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Figure  5  shows  a  plot  of 


—  — —  versus  for  —  10  for  A/  —  1,2 

9  D  rA 


Fig.  5  Probability  of  detection  versus  signal-to-noise  ratio  for 
Swerling  case  #2  target  with  N=l,2. 


2. 5  Swerling  Case  #3: 

We  see  from  figure  2,  that  the  probability  density  function  when  one  large 
scatterer  exists  with  a  number  of  small  scatterers  is,  (Nathanson,  page  86) 


p{x)  = 


X 


a 


(54) 


Substituting  the  last  line  of  equation  (54)  into  equation  (13)  gives 


4  -  J±J^] 

O  (03)  =  —— ^  e  ‘-'“W 

^  xlO-i(oy{ 

4 


(l-icoy 


^  2  imN  '' 
kK  1-1(0^ 
4 


(1  -  ico)"'"  [2  -  i(x>(2  +  Ate, )] 

1 


(1-/C0) 


N-2 


1  —  /co  1  + 


!K] 

2  ) 


“|2 


(55) 


which  agrees  with  Fehlner,  page  26,  his  equation  (A-8).  Substituting  equation  (55)  into 
equation  (14)  gives 
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pAY)  = 


2/1 1  -(1 

l  2  j 


d(0 


(56) 


“  Nx 


1  + 


Again,  we  consider  various  cases  for  N  starting  with  N=l: 


P\(X)~  f  n2  j~?  ^  ,.2  d(o 


27c|^1  +  ^ 

1_ 

27c|  1+^ 

2 

1 


1 


-iOiY 


(o),  -  i(£>y 


00 

dca-ij 


(He 


-tcoK 


2k 


fi  jc  Y 

1+— 

2j 


f 


-CD  + 


-~(co^  -  ico)' 


din 


V 


y)\ 


Ik 


t 


^e-"-''[l  +  y(l-co,)] 


.  1+—  . 

2) 


1 


X  ^ 

2k\  1  +  ^ 

2J 


2) 


1  +  -^ 

2 


(57) 


and  the  corresponding  probability  of  detection  is 
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P,  =  \p,{Y)dY 

Yo 


r  ^ 


Y 

U  J 

1+^ 

2 


dY 


(58) 


An  interesting  observation  for  this  particular  case  is  that  when  =  1,  Fehner’s 
derivation,  (his  equation  A-3)  yields  infinity  in  the  denominator,  (unity  divided  by  infinity 
equals  zero),  and  this  is  the  reason  we  suggest  using  equation  (45).  The  mathematical 
result  is  correct,  however,  it  is  not  convenient  for  numerical  implementation. 

The  probability  of  false  alarm  for  this  case  is 


(59) 


N=2: 


-  J__J_  f  g 

27C  (l  +  y  i(co^  -  iaif 


d(0 


The  probability  of  detection  is 


(60) 
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=  e  '*'•  1  + 


1  +  X 


and  the  corresponding  probability  of  false  alarm  is 


p  =e-Hi+r) 


N=3: 


PAY)  = 


j-  _ i _  r _ _ 

I' j  ^  3x^  J  i(l  -  z(o)(co^  -  my 


(1  -  ico)(co„  -  icof  1  -  m  (©,  -  mf  co^  -  m 
1  =  a(co^  -  mf  +  5(1  -  m)  +  C(1  -  icc))(co^  -  m) 
-coHA  +  C)  =  0 

-(o(2ico^A  +  iB  +  /Cco^  +  iC)  =  0 
Aco^  +5+ Cco  =  1 

a  a 

A  =  -C 


and  continuing. 
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1 


(62) 


The  probability  of  detection  becomes 


2 


and  the  probability  of  false  alarm  is 

p„=(i+r>-'-.  («) 


The  analysis  required  to  solve  for  Y  for  N— 2,  is  as  follows. 
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log.  PfA  =  log^  (l  + 


Y:  +  log,  P,X  +21og,p,,=0 _ 

1"„  =  -^log.  Pm  +  -^^Se^FA 


Figure  6  shows  a  plot  of  =  -7  versus  P^  for  P^^  =  10  ,  for 

N  =  1,  2  and  3 .  in  order  to  be  able  to  compare  our  results  with  DiFranco  and  Rubin, 
their  definition  for  normalized  signal-to-noise  power  has  to  be  related  to  our  definition  in 
equation  (7).  DiFranco  and  Rubin  define 


X 


D-R 


This  results  in  the  mean  signal  power  power  as 


A^  = 


JC  = 


Exp{xl,}  =  -A^, 
pV2  _ 


2A' 
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Fig.  6  Probability  of  detection  versus  signal-to-noise  ratio  for 
Swerling  #3  target. 

For  example,  from  figure  6 


P^=0.1,V  =  l,x  =9.0<iB. 

From  figure  11.4-1,  page  411  of  DiFranco  and  Rubin  at 


P,  =  0.1  (10%),  91^  =  12.2  -3.0  =  9.2  , 

where  we  need  to  subtract  3  dB  because  DiFranco  and  Rubin  use  peak  power  and  we  use 
average  power.  The  above  results  are  within  numerical  accuracy  in  reading  their  plot. 

From  Fehlner,  page  67,  figure  31,  for  =  0. 1,  S  /  N  =  8. 0  (9. 0  dB).  For 
N  =  2,  0-1,  from  figure  6,S  I  N  =  6.5  dB,  and  from  Fehlner,  figure  32, 

page,  iV  =2,  P^=0.l,S/  N  =  6.53  dB,  which  is  excellent  agreement. 
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2.6  Swerling  Case  #4: 

As  discussed  in  section  1.1,  the  order  of  applying  the  pdf  for  the  target  power,  and 
treating  the  sum  of  N-puises  needs  to  be  interchanged  from  the  previous  Swerling  cases. 
To  evaluate  the  characteristic  function  for  SC-4,  we  average  with  respect  to  p(x).  For 
SC-4,  pix)  is  the  pdf  for  one  dominant  scatterer  plus  noise.  For  SC-2,  p(x)  is  the  pdf 
for  all  the  returns  in  figure  2.  For  SC-4,  we  first  average  over  pix)  and  then  address  the 
N-pulse  case.  Using  the  same  pdf  for  the  target  as  in  SC-3  (Nathanson),  we  have  from 
the  last  line  in  equation  (9), 


0(0))  =  J  <E)(co|x)  p(x)dx 

0 


1 


1- 


(l-zco)of  X 


2j 
(1  -  /co) 


-I2 


2J 


(65) 


Substituting  equation  (65)  into  equation  (14)  (O^  (CO)  —  [O(C0)]  )  gives  for  the  N- 
pulse  case 


pAY)  = 


1 


1 


,2N 


I  2 


Jp - ^ 


1 


1  +  -^ 

2 


tco 


(66) 


where  the  integrand  of  equation  (66)  agrees  with  Fehlner  (1962,  page  426,  equation  A- 
13). 
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We  again  consider  cases  beginning  with  N=l: 


pxy)= 


2k  ( .  X 
1  +  ^ 

I  2 


j 


1+^ 

2 


-ICO 


=  coV““‘'  l  +  co  ^  F 


which  agrees  with  SC-3  for  N=l.  The  corresponding  probability  of  detection  is 


F^  =  \p,{Y)dY 


=  l  +  col^jF 


CO  = 


The  probability  of  false  alarm  is 


P 

^  FA 


N=2: 
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P.iY)  = 


1 


1 


2 


oo 


(1  — ico)  e 


2 


f 


\ 


d(0 


1 


1  +  ^ 
2 


-/CO 


y 


1  1  7(1-2/co-co-) 


271 


V 


1+-^ 

2 


j 


r 


1 


1+— 

2 


/CO 


=  coT^e"““" 

a 

1 


CO  = 


^  CO^X^  x,co 


J 

n 


V 


y 


l  +  _iL 

2 


The  probability  of  detection  is 


and  the  probability  of  false  alarm  is 


p„=e-Hi+i:) 


(72) 


N=3: 


P.(Y)  = 


1  1  f  . 

- —d(0 


2 


V 


^ 


1 


Y^e 


r(6)|i+^ 


V 

1+ 

r 


t  X 
1+-^ 

2 


m 


J 


-00  +  — 
I  “ 


+ 


co^ 


Y) 
lOco  20 


+ 


V 


-co^ 


Y  Y\ 
15co^  60co 


60 


(73) 


The  probability  of  detection  is 


aj 


(74) 


The  probability  of  false  alarm  is 
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Fig.  7  Probability  of  detection  versus  signal-to-noise  ratio  for 
Swerling  case  #4  target. 


2.7  Swerling  Target  plus  Compound-k  Clutter  and  Gaussian  Noise 


The  probability  density  function  for  the  sum  of  the  target  plus  clutter  signal  in 
equation  (3)  is 


(76) 


For  a  Rayleigh  target,  with  uniformly  distributed  phase,  (D,  we  know  that 


p,(x,(0) 


(77) 


and  the  pdf  for  the  sum  in  equation  (76)  of  two  Rayleigh  phasors  (i.e.,  y  is  also  Rayleigh 
distributed)  is  (Beckmaim,  page  124,  equation  4.5-4) 


0 


/ 

V 


PY 


A 


O 


m- 


T  J 


(78) 


In  section  4,  we  show  that  when  a  Rayleigh  pdf  representing  “spikey”  clutter, 
with  clutter  variance,0^,  is  modulated  by  a  gamma  pdf,  the  resulting  pdf  is  the  so-called 
compound-K  clutter  pdf.  The  principal  features  of  a  high-resolution  radar  return  from  the 
sea  consist  of  a  fast  time  variation  within  a  resolution  cell  (the  so-called  temporal 
variation  on  time  scales  on  the  order  of  several  millisec)  and  a  slow  time  variation  from 
cell  to  cell  (the  so-called  spatical  variation  on  time  scales  on  the  order  of  several 
seconds).  We  assume  the  sea  is  frozen  on  a  pulse-to-pulse  basis  (this  would  be  valid  for 
pulse  repitition  intervals  of  10  millisec  based  on  the  decorrelation  time  for  sea  scatter). 
The  fast  time  features  of  the  return  are  “spikey”  and  have  a  Rayleigh  pdf.  The  reason  the 
short  time  statistics  are  Rayleigh,  is  the  surface  slopes  on  the  sea  follow  a  Gaussian  pdf, 
and  the  resulting  scattered  field  is  Gaussing  in  the  I  and  Q  components,  or  a  Rayleigh 
eavetepe.  The  slow  time  variation,  or  the  envelope  of  the  return  is  modulated  by  the 
gamma  pdf,  because  the  swell  moves  in  and  out  of  a  particular  slant  range  cell.  As  shown 
in  section  4,  when  the  compound-K  statistics  are  included,  we  will  require  the  clutter 
variance  in  the  Rayleigh  pdf  to  satisfy 


41 


k"  +3k'  -2k  +  1 


(79) 


where  OL  defines  the  mean  value  of  the  clutter  power  and  K  defines  a  shape  parameter 
that  defines  the  clutter  variance.  In  order  to  proceed,  we  assume  in  (78)  that 


(80) 


where  the  pdf  for  the  clutter  amplitude  is  conditioned  on  the  clutter  variance  as  shown  in 
section  4.  As  will  be  seen  in  section  4,  this  conditional  form  for  the  pdf  introduces 
mathematical  complications  which  do  not  permit  writing  the  probability  of  detection  in 
closed  form,  and  to  proceed,  we  take  the  approach  that  the  clutter  variance  in  (80)  be 
treated  as  simply  an  ordinary  probability  assuming  O]  is  a  variable.  Substituting  (80) 

into  (78)  yields 


pV2 

Recall  that  X  — - 


,  and  it  is  easy  to  show  that 


X 
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We  now  substitute  (82)  in  (13)  to  obtain 


<1>„(0))  = 


(0 


(0  = 


(1-iCO)^  ‘(co^  -/©) 
1 


l  +  N  ^  +  ^ 


Now  substitute  (83)  into  (14)  to  find 


2k  i(l  -  i(oY~'  (co^  -  ico) 


d(0 


where,  in  several  steps  above  we  have  made  use  of  the  following  result 


F(a,  b)  =  \  /„  {bx)xdx 

0 


L(.bx)  =  t, 


4  J 


J.- 


r:ik'.r(k  +  l) 

1 


'  = 


2a 


2k+2 


n/sr+i). 


We  now  examine  specific  values  for  N. 


N=l: 


p,(y)=cox"^ 


(86) 


and  the  probability  of  detection  is 


P,=e 


(87) 


The  corresponding  probability  of  false  alarm  is 


P  =e 

^  FA 


(0"  =1  + // 


x  = 


l+N 


J\ 


HPfa) 


(88) 


N=2: 


O,(co)  = 


1 

_  1  ^ 

(l-co.) 

^00^  -fCD 

3 

1 

1 

(89) 


and  substituting  (89)  into  (14)  gives 


(0 


1-co 


(90) 
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and  the  probability  of  detection  is 


P  = 


© 

c 

r  1 

1 

1 

(l-Q),) 

- c 

l®c 

^ ) 

(91) 


and  the  probability  of  false  alarm  is  given  in  (88). 


N=3: 


P,iY)= 


i-e-'  +  Ye-''  +6-"'’') 


The  probability  of  detection  is 


p  ^ 
^  D 


(1-®  y 


(92) 


(93) 


and  the  probability  of  false  alarm  is  given  in  (88). 

In  general,  from  induction,  the  probability  of  detection  for  this  case  for  any  N  is 


-«»cn, 


(94) 


'-'r  r» 

Figure  8  shows  a  plot  of  — -  in  dB  versus  for  a  target-to-clutter  ratio. 


=  1  (a  “stressing”  case),  for  N=l,  2  and  3. 
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Pd 


Fig.  8  Probability  of  detection  versus  signal-to-noise  ratio  for  a 
target  in  clutter  plus  noise.  Three  different  cases  for  the 
number  of  integrated  pulses  are  considered. 

From  figure  8  we  see  for  the  case  examined,  where  the  target-to-clutter  ratio  is  10 
(a  stressing  case),  there  is  relatively  little  significant  improvement  by  integrating  more 
pulses. 

We  now  give  an  alternative  derivation  for  probability  of  the  sum  in  (76).  The  pdf 
for  target  plus  noise  can  be  written  as 


1 

2koI 


e 


and  we  wish  to  determine  the  pdf  of 


(95) 
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(96) 


p  =  ^/?+/ 
JC  =  T  e'“ 


y  = 


We  make  the  transformation  of  variables 


X  =  pcose,  y  =  psinG,  dxdy  =  pdprfG 


to  obtain 


1  23c  (pcose-T)^+(psin9f 

^  ^  pt/pc?e 


p,  Jp!i!!12,c 

__P_^  24 
27C0^  0 


2cl 


^2 


pT 

WrJ 


Again,  using  the  variable  X  — 


p^  /  2 


,  equation  (98)  becomes 


pWt)  = 


(<] 

(d  dn] 
^“Ur^  4  jy 

0 

v4> 

1  J 

Then  from  the  Theorem  of  Total  Probability 


(97) 


(98) 


(99) 
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pW=JpWy)p(Y)^Y 


p{x)  = 


J' 


vo;ao:/ 


A 


A 


7 

0 

<^n  00 

Q  1 

1 _  ■ 

e 

0 

X 

^y^/lxc^ 


^ 


ol 


\p{y)dy 


2(<Jj+a^) 


yv^c„ 


^ 


V 


yiy 


7 


in  agreement  with  equation  (81). 


(100) 
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3.  An  Alternative  Derivation  for  the  Probability  Density  Function  for  a  Target  and 
Gaussian  Noise  Fluctuating  Pulse-to-Pulse  (Swerling  Class  2,  N=l). 

As  an  alternative  to  the  derivation  in  section  1.2,  consider  the  following  analysis. 
When  the  target  varies  so  rapidly  that  its  cross  section  changes  from  pulse-to-pulse  along 
with  the  noise,  we  can  replace  equation  (1)  for  the  total  received  signal  in  the  I  and  Q 
channels  (no  clutter)  as 


Re'«=Te‘“  +  re'\ 


(101) 


and  the  probability  density  function  for  the  received  signal  in  equation  (101)  is,  assuming 
the  target  is  Rayleigh  distributed,  and  the  noise  is  a  uniformly  distributed  phasor  (with 
random  amplitude). 


p{F\r)  =  ^e 


2gI 


f— 1 


(102) 


X- 


p{y\x)  =  e-^^^^^lX24^),  y>0 


Then  applying  the  Theorem  of  Total  Probability, 


P{y)  =  J  p(y\x)p{x)dx,  (103) 


and,  assuming  the  noise  in  equation  (103)  is  Rayleigh  (Gaussian  in  I  and  Q  and  Rayleigh 
envelope) 
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.2 


p(x)dx  =  —e 


r^2\ 


2a^  ’ 


p{y)dy  = 


f  A 

'  <5-  ' 


of 


dy  = 


dy 


xdx 


(104) 


and  substituting  equation  (104)  into  equation  (103)  gives(x  =  ^,dz  =  2^^) 


^7  0 


1 

^rj 


(105) 


-y- 


(oj  +  o^) 


a  simple,  closed  form  result. 

An  alternative  derivation  from  Papoulis  (pages  194-195, 1965)  is  as  follows: 


=Z^  + 


p{x,r)  = 


2(0, M 


il’ 


27c(Gf  +  a„^) 


(106) 


and  assume  the  target  has  a  mean  value  of 
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Exp{'t]  =  'i^, 


(107) 


so  the  pdf  in  equation  (106)  becomes  the  Rice  pdf  as 


2n[G]+(5l)i 


(zcosi3~Ti)"+(zsinT3)~ 

2(0? +aj) 


(of  +  oD* 


2(0? +oi) 


(108) 


We  are  again  interested  in  obtaining  the  probability  density  of  the  square  of  the  envelope 
of  the  total  received  signal.  Thus  define 


^~2(of  +  c=)’^“2(af+a.") 

A%  Zt/Z 
^  =  T~2 - JY’ 


(109) 


Now,  we  need  to  integrate  the  last  line  of  equation  (109)  over  all  possible  values  of 
?{x,)as 


P(l)  =  i  e%(2M)dt; 

0  0 

^  =  x^,dt^  =  2xdx 

p{k)  =  2e-^]e-^'l„(2^x^)xdx  =  e-^=e~^ 


(110) 


which  agrees  with  equation  (104),  with  the  following  substitutions 
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2(a^+a„-)’^  2gJ 


p{z)dz  =  p{^)d^ 

Piz)  =  p(^)^- 


2(0; +4)  ’ 


dz  (a; +09 


y  =  —,dy  =  zdz 


(111) 


From  equation  (105)  we  can  derive  the  probability  of  detection,  and 
probability  of  false  alarm,  .  We  find 


yo 

00 

P.^^le-^dy^e'^-  (112) 

>'c=-log,(P^J 


In  figure  9  we  plot  the  probability  of  detection  for  three  false  alarm  numbers,  or 
probability  of  false  alarms  determined  from 


FA 


0.693147 


(113) 


where  n'  is  the  false  alarm  number.  The  curves  in  figure  9  match  those  given  by  Fehlner 
(1962). 
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4.  Probability  Density  Function  for  Broad-Band  Clutter  Constant  Pulse-to-Pulse 
but  Varying  Scan-to-Scan  (No  Target,  Probability  of  False  Alarm) 


From  the  free-space  radar  range  equation,  the  received  clutter-to-noise  power  is 


N  {4KfR‘(kTj 

and  the  clutter  variance,  ,  in  the  Rayleigh  pdf. 


(114) 


(115) 


is  given  by 


al  =  Exp{o}  =  .  (116) 


Ward  et  al.(1990)  and  Baker  (1988),  treat  clutter  by  assuming  the  clutter  is 
Rayleigh  distributed  on  time  scales  of  millsec  as  in  equation  (80),  and  assume  the 
variance  of  the  clutter  in  equation  (1 15)  is  modulated  with  a  Gamma  probability  density 
function,  on  time  scales  of  seconds,  as 


r(K+i) 


T1 


(117) 


where  the  random  variable  Tj  (the  clutter  power)  in  equation  (1 17)  equals  the  clutter 
variance  according  to 


Tl  =  2a\ 


(118) 
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and  a  is  the  mean  value  of  the  clutter  power  and  K  is  a  shape  parameter  that  defines  the 
normalized  clutter  variance. 

The  probability  density  function  for  the  clutter  amplitude,  y,  is  from  the  Theorem 
of  Total  Probability 


OQ 

/’(T)=Jp.(Yl'n)p(Ti)rfn 


2a’'^’y  f  ^  ‘ 

— - —\'x\  e'^—e  Vrj 

r(K  +  l)o  T| 


(119) 


The  integral  in  equation  (119)  is  evaluated  as  follows: 


arj +  —  =  2  cosh  (0 

ari"  -  2ti  Vay  cosh  co  +  y '  =  0 
2riy_ 


Va 


coshcD  +  y  =0 


T|  =  cosh  CO  ± ..  —  cosh^  CO  - — 

va  V  a  a 

y 

= -7=(cosh  CO  ±  sinh  co) 

Va 


->oo, 


— CX3* 


->oo,(+) 


•^o,(+) 


dT[  =  Tj  Jco 


r  Y 


0  Ti 


4B.) 


J' 


-2*N/aj:coshtD 


(cosh  CO  +  sinh  co)"  c/co 


(120) 
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The  last  integral  is  evaluated  for  the  following  cases: 
K  =  l: 


Y  f  -ij^ircoshd) 


2x 


(cosh  (0  +  sinh  co)  Joo  =  -7=/  e 

Vex  0 


-2  Vox  cosh  0) 


coshco<ia) 


2x 

Va 


KX24ax) 


k  =  2: 

^jlY 


-“2-yarcosha) 


(cosh  CO  +  sinh  co)^  <ico  =  2 


r 


-4=\  fe'^’^^'“''“(cosh^co  +  sinh^0)t/© 

VoCy  0 


=  2 


r 


-^1  L-^'^^“''’“cosh2cot/co 
va;  0 


4:,(2-x/ay) 


k  =  3: 


Y  I  f  -2A/aTCOsh(D 

i 


Z'  y  \ 


(cosh(0  +  sinhco)^Jco  =  2  — ^  k  {2^X) 

\^aj 


(121) 


and,  in  general, 


0 


T1 


(122) 
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and  from  (119)  and  (122), 


4a''*' y  y 

r(K  +  l)l^Va 


KS2^y) 


4a^'y’'*' 

r(K+i) 


(120) 


Figure  10  shows  a  plot  of  the  pdf  in  equation  (120)  together  with  a  Rayleigh  pdf 
for  comparison.  If  the  form  of  the  compound-K  pdf  given  in  equation  (1 19)  is  substituted 
in  the  right-hand-side  of  equation  (78)  for  the  resulting  integrand  containing  the 

product  of  two  Bessel  functions  does  not  lend  itself  to  a  closed  form  representation. 

From  figure  10,  we  see  the  tails  of  the  compound-k  distribution  extend  further  out  than 
the  corresponding  tails  for  the  Rayleigh  pdf.  This  probably  accounts  for  the  fact  that  the 
compound-k  pdf  predicts  more  “spikey”  clutter  than  the  Rayleigh  pdf.  This  is  also  the 
reason  in  the  design  of  a  broad-band  radar  where  the  slant  range  cell  is  small,  compared 
to  that  for  a  narrow  band  radar,  that  even  though  the  slant  range  cell  size  is  small,  the 
clutter  may  appear  more  spikey  than  the  corresponding  return  from  a  narrow  band  radar. 

Figure  10  also  shows  the  Rayleigh  pdf  over  the  same  range  of  the  independent 
variable.  The  mean  value  of  the  clutter  parameter,  Ct ,  is  adjusted  so  that  both  pdfs  have 
equal  variances.  We  will  see  in  equation  (127)  this  is  accomplished  when  K  =  1,  if 

V3 
a  =  — . 

a 

c 
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Fig.  10  A  comparison  of  the  compound-K  and  Rayleigh  pdf  s. 


We  now  show  a  derivation  for  an  approximate  form  for  the  compound-k  pdf. 
Using  the  following  asymptotic  representation  for  the  modified  Bessel  function  in 
equation  (123) 


(124) 


equation  (123)  can  be  written  as 


2jna2j__ 
r(K  +  l) 


(125) 
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and  making  the  replacement  in  (125) 


1 

K  +  — =  V, 

2 


(126) 


gives  a  form  for  p{y)  identical  to  the  gamma  pdf  (Papoulis,  1965,  p.  103-104)  as 


(127) 


with 


(128) 


Thus,  the  probability  density  function  for  the  clutter  voltage  (no  target)  can  be 
approximated  as 


r(K+i) 


(129) 


and  we  are  back  to  a  simple  gamma  pdf.  The  pdf  in  equation  (129)  represents  a  much 
simpler  expression  for  use  in  integrals  involving  the  Theorem  of  Total  Probability.  In 
figure  1 1 ,  we  show  a  plot  of  the  approximate  compound-k  (gamma)  pdf  togther  with  the 
exact  pdf. 
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Fig.  1 1  Comparison  of  the  approximate  compound-k  (equation  125)  and 
the  exact  compound-k  pdf  (equation  1 19). 

We  have  shown  in  figure  1 1  that  using  the  approximate  form  for  the  compound-k 
pdf  differs  insignificantly  from  the  exact  expression,  and  therefore  this  simpler  form 
should  be  considered  in  derivations  requiring  the  compound-k  pdf. 


The  first  two  moments  of  Tj  are  given  by 


£x/7{y}  =  Y=  ^ 


Ka 


Ea:p{(y-y)'}  =  <=-^— ify-^1  YV-"<iY  . 


r(K+i)ov  Kay 
(k'^  +3k^  -2k  +  1) 
(koT 


(130) 
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The  clutter  variance,  from  equation  (130)  is  defined  as 


(k^+3k^-2k  +  1) 

(aK)"* 


(131) 


The  clutter  variance,G^,  has  been  estimated  to  be  0.7  for  low-sea-state  clutter,  1.1 

for  ground  clutter  and  as  high  as  1.6  for  high-sea-state  clutter.  Table  1  shows  the 
corresponding  values  for  OC,  in  the  compound-K  pdf 

Table  1.  Clutter  variance  and  mean  clutter  power  for  various 
environments. 


environment 

a: 

a" 

low-sea-states 

0.7 

4.29 

ground 

1.1 

2.73 

high-sea-state 

1.6 

1.875 

The  moments  for  the  gamma  distribution  can  be  expressed  in  terms  of  the 
characteristic  function  as  (Papoulis,  1965) 


1  mi 

O^(co)  3©  ’ 


(132) 


and,  if  the  moments  in  equation  (133)  are  normalized  to  the  variance.  Ward  has  shown 
that  all  moments  greater  than  i  =  2  are  zero  (Ward,  1990).  This  reduces  to  a  Gaussian 
distribution  using  the  Central  Limit  Theorem  argument. 
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6.  Computer  Code 


c 

C  CALCULATE  THE  PROBALITY  OF  DETECTION  FOR  SC-0 
C 

COMMON/SCO/B 

OPEN(UNIT=10,FILE=’C:\tecplot\scOLdat’) 

OPEN(UNIT=l  l,FILE=’C:\tecplot\sc02.dat’) 

Y0=8*ALOG(10.) 

DO  1 1=1,40 
n=i 

B=n 

Z  =  2.*SQRT(B*Y0) 

CALL  TRAP1(0.,Z,SUM,15) 

PDl  =  1.  -  EXP(-B)*SUM 
WRITE(10,*)  PDLIO  *ALOG10(B) 

WRITE(*,*)  PD1,B 

1  CONTINUE 
A=  Le-08 

YO  =  (-0.5*ALOG(A)  + 

X  SQRT(0.25*ALOG(A)**2  -  2.*ALOG(A))) 

DO2J=l,40 
FJ  =  J 
B  =  FJ 

Z2  =  2  *SQRT(B*Y0) 

CALL  TRAP2(0.,Z2,SUM2,15) 

PD2  =  1.  -  EXP(-2.*B)*SUM2 
WRITE(I1,*)  PD2,I0.*ALOG10(B) 

WRITE(*,*)  PD2,B 

2  CONTINUE 
STOP 
END 


SUBROUTINE  TRAP1(A,B,SUM,N) 
EXTERNAL  FSCO 
NN=N+ 1 
SUM  =  0. 

FN  =  N 
D  =  (B-A)/FN 
X  =  A 

DO  10I=1,NN 
J=  (I-1)*(NN-I) 

C=  1. 

IF(J)  15,5,6 
5C  =  0.5 

6  SUM  =  SUM  +  C*FSC0{X) 

10X  =  X  +  D 
SUM  =  SUM*D 
RETURN 

15  WRITEC*  *)  ’INCORRECT  INDEX’ 
STOP 
END 
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non  n  n  n 


SUBROUTINE  TRAP2(A,B,SUM,N) 
EXTERNAL  FSC02 
NN=N+ I 
SUM  =  0. 

FN  =  N 
D  =  (B-A)/FN 
X  =  A 

DO  10I=1,NN 
J=  (I-1)*(NN-I) 

C=  1. 

IF(J)  15,5,6 

5  C  =  0.5 

6  SUM  =  SUM  +  C*FSC02(X) 
10X=:X  +  D 

SUM  =  SUM*D 
RETURN 

15  WRITE(*,*)  ’INCORRECT  INDEX’ 
STOP 
END 


A  FUNCTION  USED  IN  SC-0 

FUNCTION  FSCO(X) 

EXTERNAL  BESSIO 
COMMON/SCO/B 

FSCO  =  BESSI0(X)*X*EXP(-X**2/(4.*B)) 

RETURN 

END 


A  FUNCTION  USED  IN  SC-0 

FUNCTION  FSC02(X) 

EXTERNAL  BESSIl 
COMMON/SCO/B 

FSC02  =  BESSI1(SQRT(2.)*X)*(X**2/(4.*SQRT(2.)*B**2))* 
X  EXP(-X**2/(4.*B)) 

RETURN 

END 


non  non 


PROBABILITY  OF  DETECTION  FOR  SC-1 

OPEN(UNIT=10,FILE=’C:\tecplot\scl.dat’) 

Y0=  10  *ALOG(10) 

DO2K=l,10 

N=K 

PWR  =  N  -  1 
DO  1 1=  1,100 
XA  =  P.2 

PD  =  ((1.  +  (1./(N*XA)))**PWR)*EXP(-Y0/(1.  +  N*XA)) 
WRITE(*,*)  PD,10.*ALOG10(XA) 

WRITE(10,*)  PD,10.*ALOG10(XA) 

1  CONTINUE 

2  CONTINUE 
RETURN 
END 


PROBABILITY  OF  DETECTION  FOR  SC-2 

EXTERNAL  FACTRL 
OPEN(UNIT=10,HLE=’C:\tecplot\sc2.dat’) 

Y0  =  8*ALOG(10) 

DO  1 1=1,100 
XA  =  1*0.2 

PD  =  EXP(-Y0/(1.  +  XA)) 

WRITE(*,*)  PD,10.*ALOG10(XA) 
WRITECIO,*)  PD,10.*ALOG10(XA) 

1  CONTINUE 
RETURN 
END 


non 


PROBABILITY  OF  DETECTION  FOR  SC-3 

EXTERNAL  FACTRL 
OPEN(UNIT=  10,FILE=’C:\tecplot\sc3.dat’ ) 

Y0  =  8*ALOG(10.) 

WRITE(*,*)  YO 
DO  1  1=1,200 
XA  =  I 

OMEGA  1  =  l./(L  +  (XA/2.)) 

PDl  =EXP(-YO*OMEGA1)*(L  +  (XA/2.)*YO*OMEGA1**2) 
WRITEdO,*)  PD1,10.*ALOG10(XA) 

C  WRITE(*,*)PD1,10.*ALOG10(XA) 

1  CONTINUE 
A=  Le-08 

Y0  =  (-0.5*ALOG(A)  + 

X  SQRT(0.25*ALOG(A)**2  -  2.*ALOG(A))) 

WRITE(*,*)  YO 
DO  2 1=1,200 
XA  =  I 

OMEGA2  =  l./(l.  +  XA) 

PD2  =  EXP(-YO*OMEGA2)*(1.  +  YO*OMEGA2) 
WRITEdO,*)  PD2,10.*ALOG10(XA) 

C  WRITE!*,*)  PD2,10.*ALOG10(XA) 

2  CONTINUE 
DO  3  1=1,200 
XA  =  I 

OMEGA3  =  l./(l.  +  (3.*XA/2.)) 

PD3  =  EXP(-Y0*OMEGA3)*(l.  +  d./(3.*XA/2.))*(l.  +  YO)) 
WRITEdO,*)  PD3,10  *ALOG10(XA) 

C  WRITE!*,*)  PD3,10.*ALOG10!XA) 

3  CONTINUE 
RETURN 
END 


non 


PROBABILITY  OF  DETECTION  FOR  SC-4 


OPEN(UNrr=10,FILE=’C:\tecplot\sc4.dat’) 

Y0  =  8*ALOG(10.) 

WRITE(*,*)  YO 
N=  1 

DO  1 1=1,200 
XA  =  I 

OMEGA  =  17(1.  +  N*(XA/2.)) 

PDl  =  EXP(-YO*OMEGA)*(1.  +  (XA/2.)*OMEGA**2*YO) 
WRITEdO,*)  PD1.10.*ALOG10(XA) 

C  WRrrE(*,*)PDl,10.*ALOG10(XA) 

1  CONTINUE 
A=  l.e-08 

Y0  =  (-0.5*ALOG(A)-h 

X  SQRT(0.25*ALOG(A)**2  -  2.*ALOG(A))) 

WRITE(*,*)  YO 
DO  2 1=1,200 
XA  =  I 

OMEGA  =  1-/(1.  -H  0.5*XA) 

PD2  =  OMEGA**2*EXP(-YO*OMEGA)*((1.  +  0.5*XA)**2  + 
X  OMEGA*YO*(  1 .  0.5*XA)**2  + 

X  0.5*XA*OMEGA**2*Y0**2*(l.  +  0.25*XA)  + 

X  (176.)*(0.5*XA)**2*(Y0*OMEGA)**3) 

WRnE(10,*)  PD2,10.*ALOG10(XA) 

C  WRrrE(*,*)PD2,10.*ALOG10(XA) 

2  CONTINUE 
PFA=  l.e-08 

CALL  NEWTON(PFA,YO) 

WRnE(*,*)  YO 
DO  3  1=1,200 
XA  =  I 

OMEGA  =  17(1.  -H  0.5*XA) 

PD3  =  (0.5*XA)*(OMEGA**6)*EXP(-Y0*OMEGA)/120.* 

X  (Y0**5*(l.  -  2.*OMEGA  +  OMEGA**2) 

X  +(57OMEGA)*Y0**4*(l .  +  OMEGA  -  2.*OMEGA**2) 

X  +(20./OMEGA**2)*(1.  -i-  OMEGA  -t-  OMEGA**2)*( 

X  Y0**3  +  (37OMEGA)*Y0**2  +  (6./OMEGA**2)*Y0  + 

X  (6./OMEGA**3))) 

WRITE(10,*)  PD3,10.*ALOG10(XA) 

C  WRITE(*,*)PD3,10.*ALOG10(XA) 

3  CONTINUE 
RETURN 
END 
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n  n  n 


FIND  THE  THRESHOLD  USING  NEWTON’S  METHOD  FOR  SC4 

SUBROUTINE  NEWTON(PFA,YO) 

YOI  =  -ALOG(PFA) 

R1  =  (1.  +  YOI  +  0.5*Y0I**2) 

R2  =  EXP(-YOI) 

R3  =  O.5*Y0I**2 

YO  =  YOI  -  (PFA/(R2*R3)  -  R1/R3) 

RETURN 

END 
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non 


PROBABILITY  OF  DETECTION  FOR  SC-5 

OPEN(UNrr=10,FILE=’C;\tecplot\sc5.dat’) 

A  =  lO.e-08 
N=1 

XB  =  10. 

YO  =  -(1.  +  N*XB)*ALOG(A) 

WRITE(*,*)  YO 
DO  1 1=1,200 
XA  =  20.*I 

OMEGACl  =  17(1.  +  (XB  +  XA)) 

PDl  =EXP(-OMEGAC1*YO) 

WRITE(10,*)  PD1,10.*ALOG10(XA) 

WRITE(*,*)  PD1,10.*ALOG10(XA) 

1  CONTINUE 
N=2 

YO  =  -(!.  +  N*XB)*ALOG(A) 

DO  2 1=1,200 
XA  =  20.*I 

OMEGAC2  =  l./(l.  +  2.*(XB  +  XA)) 

PD2  =  (l./(l.  -  OMEGAC2))*EXP(-OMEGAC2*YO) 
WRITE(10,*)  PD2,10.*ALOG10(XA) 

WRITE(*,*)  PD2,10.*ALOG10(XA) 

2  CONTINUE 
N=3 

YO  =  -(1.  +  N*XB)*ALOG(A) 

DO  3  1=1,200 
XA  =  20.*I 

OMEGAC3  =  l./(l.  +  3.*(XB  +  XA)) 

PD3  =  (!./((!.  -  OMEGAC3)*»2))*EXP(-OMEGAC3*YO) 
WRITE(10,*)  PD3,10.*ALOG10(XA) 

C  WRITE(*,*)PD3,10.*ALOG10(XA) 

3  CONTINUE 
RETURN 
END 
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c 

C  PDFFORTHECOMPOUND-K 
C 

EXTERNAL  FUNCTION  BESSKl 
OPEN(UNrT=10,FILE=’C:\tecplot\compound.dat’) 

OPEN(UNrr=  1 1  ,FILE=’C:\tecplot\ray  Ieigh.dat’ ) 

OPEN(UNIT=  12,FILE=’C:\tecplot\approxim.dat’ ) 

PI  =  3.1415926536 
ALPHA  =1. 

KAPPA  =  1. 

SIGMAC  =  SQRT(3.)/ALPHA 
PWRl  =  (1/2)  +  KAPPA 
PWR2  =  2 
DO  1 1=1,100 
n  =  i 

x  =  n*o.i5 

COMPK  =  4.*(ALPHA**PWR1)*(X**PWR2)*BESSK1(2.*SQRT(ALPHA)*X) 
WRITE(10.*)  X,COMPK 

1  CONTINUE 
DO2J=l,I00 
FJ  =  J 

XX  =  FJ*0.1 

RAYLEIGH  =  (XX/SIGMAC)*EXP(-XX**2/(2.*SIGMAC)) 

WRITE(1 1,*)  XX,RAYLEIGH 

2  CONTINUE 

PWRl  =  KAPPA  +  0.5 
DO  3  K  =1,100 
FK  =  K 
XG  =  0.1*FK 

APPROX  =  2.*SQRT(PI)*XG**PWR1*EXP(-2.*XG) 

WRrTE(12,*)  XG, APPROX 
WRITE/*,*)  XG,APPROX 

3  CONTINUE 
STOP 
END 


u  u  u 


CALCULATE  THE  RICE  PDF 
EXTERNAL  FUNCTION  BESSIO 

OPEN(UNrr=  1 0,FILE=’C:\fortran\probilty\code\ricepdf.out’ ) 
OPEN(UNrr=l  1  ,FILE=’C:\fortran\probilty\code\gausspdf.out’) 
OPEN(UNIT=12,FILE=’C:\fortran\probilty\code\swerlpdf.out’) 
X=0.01 
XRCS=  1. 

SIGMA  =  0.4 
PI  =  3.1415926536 
DO  1 1=1,101 
H  =  I 

Y  =  (FI-l.)*.l 
Z  =  X*Y 

R1  =  BESSI0(2.*SQRT(Z)) 

R2  =  exp(-(X+Y))*R1 

R3  =  (l./(SIGMA*SQRT(2.*PI)))*EXP(-(Y**2/(2.*SIGMA**2))) 
R4  =  (l./XRCS)*EXP(-Y/XRCS) 

WRITEdO,*)  Y,R2 
WRITE(11,*)  Y,R3 
WRITE(12,*)  Y,R4 
1  CONTINUE 
STOP 
END 


FUNCTION  bessiO(x) 

REAL  bessiO,x 
REAL  ax 

DOUBLE  PRECISION  pl,p2,p3,p4,p5,p6,p7,ql,q2,q3,q4,q5.q6,q7,q8,q9,y 
SAVE  p  I,p2,p3,p4,p5,p6,p7,q  1  ,q2,q3,q4,q5,q6,q7  ,q8,q9 
DATApl,p2,p3,p4,p5,p6,p7/1.0d0,3.5156229d0,3.0899424d0, 
*1.2067492d0,0.2659732d0,0.360768d-l,0.45813d-2/ 

DATA  q  1  ,q2,q3,q4,q5,q6,q7,q8,q9/0.39894228d0,0. 1328592d- 1 , 
*0.225319d-2,-0.157565d-2,0.916281d-2,-0.2057706d-l,0.2635537d-l, 

*-0. 1647633d-l,0.392377d-2/ 
if  (abs(x).lt.3.75)  then 
y=(x/3.75)**2 

bessi0=p  1  +y*(p2+y*(p3+y*(p4+y*(p5+y*(p6+y*p7))))) 
else 

ax=abs(x) 

y=3.75/ax 

bessi0=(exp(ax)/sqrt(ax))*(ql+y*(q2+y*(q3+y*(q4+y*(q5+y*(q6+y* 

*(q7+y*(q8+y*q9)))))))) 

endif 

return 

END 

C  (C)  Copr.  1986-92  Numerical  Recipes  Software  L91y. 
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non 


CALCULATE  PROBABILITY  OF  FALSE  ALARM 

EXTERNAL  F 
EXTERNAL  G 
EXTERNAL  FACTLN 
COMMON/POWER/  N 

OPEN(UNTr=10,FILE=’C;\fortran\probilty\data\data.in’) 
OPEN(UNIT=ll,FILE=’C:\foitran\probilty\code\falsalrm.out’) 
READ(10,*)  CLTRN0IS,AK,N 
DO  1 1=1,101 
n  =  i 

TON  =  CFI-1.)*0.05 
Rl  =  EXP(-TON) 

R2  =  CLTRNOIS*(AK**2)*(AK+l.)*(AK+2.)/(AK**4  +  3.*(AK**3) 
X  -2.*AK+1.) 

CALL  TRAP(F,TON,20.,ANS1,10) 

CALL  TRAP(G,TON,20.,ANS2,10) 

R3  =  FACTRL(N) 

R4  =  ANS1/R3 

R5  =  (N**2)*R2*ANS1/R3 

R6  =  -N*R2*ANS2/R3 

PFA  =  R4  +  R5  +  R6 

IF(PFA  .GE.  0.)  THEN 

PFA  =  PFA 

ELSE 

PFA  =  0. 

ENDIF 

WRITE(*,*)  PFA 
WRnE(ll,*)TON,PFA 
1  CONTINUE 
STOP 
END 


FUNCTION  F(X) 
COMMON/POWER/  N 
F  =  (X**(N-1))*EXP(-X) 
RETURN 
END 


FUNCTION  G(X) 
COMMON/POWER/  N 
G  =  (X**N)*EXP(-X) 
RETURN 
END 
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FUNCTION  factrl(n) 

INTEGER n 
REAL  factrl 
CU  USES  gammln 
INTEGER  j.ntop 
REAL  a(33),gammln 
SAVE  ntop,a 
DATA  ntop,a(l)/0,U 
if  (n.lt.O)  then 

pause  ’negative  factorial  in  factrl’ 
else  if  (n.Ie.ntop)  then 
factrl=a(n+l) 
else  if  (n.le.32)  then 
do  1 1  j=ntop+l,n 
a(j+l)=j*a(j) 

1 1  continue 
ntop=n 
factrl=a(n+l) 
else 

factrl=exp(gainmln(n+ 1 .)) 
endif 
return 
END 

C  (C)  Copr.  1986-92  Numerical  Recipes  Software  i_91y. 


FUNCTION  gammln(xx) 

REAL  gammln, XX 
INTEGER  j 

DOUBLE  PRECISION  ser,stp,tmp,x,y,cof(6) 

SAVE  cof,stp 

DATA  cof,stp/76. 1 8009172947 146d0,-86.5053203294 1 677d0, 

*24.01409824083091d0,- 1 .23 1739572450 155d0,.  1 2086509738661 79d-2, 

*-.5395239384953d-5,2.5066282746310005d0/ 

x=xx 

y=x 

tmp=x+5.5d0 

tmp=(x+0.5d0)*log(tmp)-tmp 
ser=l  .000000000190015d0 
do  11  j=l,6 
y=y+l.d0 
ser=ser+cof(j)/y 
1 1  continue 

gammln=tmp+log(stp*ser/x) 

return 

END 

C  (C)  Copr.  1986-92  Numerical  Recipes  Software  i_91y. 


SUBROUTINE  trap(F,a,b,s,n) 

INTEGER  n 
REAL  a,b 
REAL  S,F,SUM 
EXTERNAL  F 
INTEGER  itj 
REAL  del,tnni,x 
if  (n.eq.l)  then 
s=0.5*(b-a)*(F(a)+F(b)) 
else 

it=2**(n-2) 

tnm=it 

del  =  (b-a)/tnm 
X  =  a  +  0.5*del 
sum  =  0. 
do  11  j=l,it 
sum  =  sum  +  F(x) 
x=x+del 
1 1  continue 

s=0.5  *  (s+(b-a)  *sum/tnm) 
endif 
return 
END 

C  (C)  Copr.  1986-92  Numerical  Recipes  Software  L91y. 


FUNCTION  besskl(x) 

REAL  besskl,x 
CU  USESbessil 
REAL  bessil 

DOUBLE  PRECISION  pl,p2,p3,p4,p5,p6,p7,ql,q2,q3,q4,q5,q6,q7,y 

SAVE  p  1  ,p2,p3 ,p4,p5,p6,p7,q  1  ,q2,q3,q4,q5,q6,q7 

DATA  p  1  ,p2,p3,p4,p5,p6,p7/l  .0d0,0. 1 5443 144d0,-0.67278579d0, 

*-0. 1 8 1 56897d0,-0. 1 9 1 9402d- 1  ,-0. 1 1 0404d-2,-0.4686d-4/ 

DATA  q  1  ,q2,q3,q4,q5,q6,q7/l  .2533 14 14d0,0.234986 1 9d0,-0.3655620d- 1 , 

*0.1504268d-l,-0.780353d-2,0.325614d-2,-0.68245d-3/ 

if  (x.le.2.0)  then 

y=x*x/4.0 

besskl=(log(x/2.0)*bessil(x))+(L0/x)*(pl+y*(p2+y*(p3+y*(p4+y* 

*(p5+y*(p6+y*p7)))))) 

else 

y=2.0/x 

besskl=(exp(-x)/sqrt(x))*(ql+y*(q2+y*(q3+y*(q4+y*(q5+y*(q6*fy* 

*q7)))))) 

endif 

return 

END 

C  (C>CopT.  1986-92  Numerical  Recipes  Software  L91y. 


FUNCTION  bessiO(x) 


REAL  bessiO,x 
REAL  ax 

DOLOBLE  PRECISION  pl,p2,p3,p4,p5,p6,p7,ql,q2,q3,q4,q5,q6,q7,q8,q9,y 
SAVE  p  1  ,p2,p3,p4,p5,p6,p7,q  1  ,q2,q3,q4,q5,q6,q7  ,q8,q9 
DATApl,p2,p3,p4,p5,p6,p7/1.0d0,3.5156229d0,3.0899424d0, 

*  1 .2067492d0,0.2659732d0,0.360768d- 1 ,0.458 1 3d-2/ 
DATAql,q2,q3,q4,q5,q6,q7,q8,q9/0.39894228d0,0.1328592d-l, 

*0.2253 1 9d-2,-0. 1 57565d-2,0.9 1 628 1 d-2,-0.2057706d- 1 ,0.2635537d- 1 , 

*-0. 1 647633d- 1 ,0.392377d-2/ 
if  (abs(x).lt.3.75)  then 
y=(x/3.75)**2 

bessi0=pl-t-y*(p2+y*(p3+y*(p4+y*(p5+y*(p6+y*p7))))) 

else 

ax=abs(x) 

y=3.75/ax 

bessi0=(exp(ax)/sqrt(ax))*(ql+y*(q2+y*(q3+y*(q4+y*(q5+y*(q6+y* 

*(q7-hy*(q8+y*q9)))))))) 

endif 

return 

END 

C  (C)  Copr.  1986-92  Numerical  Recipes  Software  i_91y. 


FUNCTION  bessil(x) 

REAL  bessil,x 
REAL  ax 

DOUBLE  PRECISION  pl,p2,p3,p4,p5,p6,p7,ql,q2,q3,q4,q5,q6,q7,q8,q9,y 
SAVE  p  1  ,p2,p3,p4,p5,p6,p7,q  1  ,q2,q3,q4,q5,q6,q7  ,q8,q9 
DATA  p  1  ,p2,p3,p4,p5,p6,p7/0.5d0,0.87890594d0,0.5 1498869d0, 
*0.15084934dO,0.2658733d-l,0.301532d-2,0.32411d-3/ 

DATA  q  I,q2,q3,q4,q5,q6,q7,q8,q9/0.39894228d0,-0.3988024d-l , 
*-0.362018d-2,0.163801d-2,-0.1031555d-l,0.2282967d-l,-0.2895312d-l, 
*0. 1 787654d- 1  ,-0.420059d-2/ 
if(abs(x).lt.3.75)then 
y=(x/3.75)**2 

bessil=x*(pl+y*(p2+y*(p3+y*(p4+y*(p5+y*(p6+y*p7)))))) 

else 

ax=abs(x) 

y=3.75/ax 

bessi  l=(exp(ax)/sqrt(ax))*(q  1  +y*(q2+y*(q3+y  *(q4+y*(q5+y*(q6+y* 
*(q7+y*(q8+y*q9)))))))) 
if(x.lt.0.)bessi  1  =-bessi  1 
endif 
return 
END 

C  (C)  Copr.  1986-92  Numerical  Recipes  Software  i_91y. 
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